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We investigate spectral properties ol a lermion coupled with a massive gauge boson with a mass 
m at finite temperature (T) in the perturbation theory. The massive gauge boson is introduced 
as a U{1) gauge boson in the Stueckelberg formalism with a gauge parameter a. We find that the 
fermion spectral function has a three-peak structure for T m irrespective of the choice of the 
gauge parameter, while it tends to have one faint peak at the origin and two peaks corresponding 
to the normal fermion and anti-plasmino excitations familiar in QED in the hard thermal loop 
approximation for T ^ m. We show that our formahsm successfully describe the fermion spectral 
function in the whole T region with the correct high-T limit except for the faint peak at the origin, 
although some care is needed for choice of the gauge parameter for T ^ m. We clarify that for 
T ~ m, the fermion pole is almost independent of the gauge parameter in the one-loop order, while 
for T 3> m, the one- loop analysis is valid only for a <^ 1/g where g is the fermion-boson coupling 
constant, implying that the one- loop analysis can not be valid for large gauge parameters as in the 
unitary gauge. 

PACS numbers: ll.lO.Wx, 12.38.Mh 



I. INTRODUCTION 

It is well known that for extremely high temperature (T) where the hard thermal loop (HTL) approximation in 
QED and QCD [1-4] is valid, a fermion (quark) coupled with thermally excited gauge fields (gluons) make collective 
excitations, i.e., the normal fermion (particle) and the anti-plasmino excitation with distinct peaks in the fermion 
spectral function [4]; this feature obtained in the HTL approximation is also known to be gauge invariant in the 
sense that the fermion self-energy at one- loop order does not depend on gauge [2]. As for lower T region, a possible 
change in the spectral properties of the quark in association with chiral transition in QCD was investigated [5], 
using the Nambu-Jona-Lasinio model [6], and it is shown that the coupling with the chiral soft modes [7] make the 
quark spectral function have distinct three peaks near but above the critical temperature of chiral transition. The 
appearance of such a novel spectral function at T ^ m was later confirmed [8] for a massless fermion coupled with an 
elementary massive boson with a mass m, irrespective of the type of the massive boson. The mechanism for realizing 
the three-peak structure in the spectral function was also elucidated [8] in terms of the Landau damping owing to the 
collisions of the fermion with thermally excited bosons^. 

Then one may naturally ask a question if the fermion spectral function at those lower T would smoothly connect 
with that at extremely high T, i.e. the HTL result in QED/QCD: If it is not the case, it means that we do not have 
a unified understanding of the fermion spectral properties in the whole T region. Partly to answer this question, we 
investigate spectral properties of a fermion coupled with a massive vector boson introduced as a U{1) gauge boson 
in the (generalized) Stueckelberg formalism with a gauge parameter a [10, 11], and carefully examine their possible 
gauge dependence at T 7^ 0, at the one-loop order as in [8]. Here the spectral properties include the number of the 
fermion poles, the pole position in the complex energy plane and the spectral function in the momentum-energy plane. 
We are also interested in how the quasi-particle nature of the fermion is realized or destroyed by the coupling with a 
massive boson at finite T. 

We find that the present formalism gives a valid description of the fermion coupled with a massive vector boson 
for the whole temperature (T) region at one-loop order in a unified way; thereby we reveal the characteristics of 
the fermion spectral properties depending on the distinct T regions, i.e. (I) T <^ m, (II) T ^ m and (HI) T ^ m. 
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^ This feature that the three-peak structure arises at T ~ m is not altered even for a massive fermion with a mass as long as m j is 
not too large compared with m [9]. 
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Especially, we shall show that the fermion spectral function certainly tends to have a three-peak structure for T ^ m 
in the small momentum region with supports in the positive, zero and negative energy regions. 

The investigation of the possible gauge dependence turns out to be involved owing to the appearance of a novel mass 
scale -y/aTO, inherent in the present formalism, as well as the boson mass m and temperature T. One should remark 
here that the Proca formalism adopted in [8] is not adequate for this purpose, because this formalism corresponds to a 
special gauge with a — > oo (unitary gauge), and does not lead to the proper high-T limit, or m/T 0, which should 
be the HTL approximation in QED at one-loop level [4]. This is the reason why we have adopted the Stueckelberg 
formalism to describe the massive vector boson. We remark that although the pole position is gauge-independent in 
the exact calculation [12], a gauge- dependence of the fermion pole may appear in the perturbation theory at finite 
T in general. Since the Proca formalism corresponding to the limit a — >■ oo leads to a wrong high-T limit, there 
should exist an adequate gauge-parameter region in which the results in the perturbation theory hardly show gauge 
dependence: Indeed, we show that this is the case in the present work. 

This paper is organized as follows. In Sec. II, we formulate the U (1) gauge theory in which the gauge boson acquires 
finite mass. We perform a calculation of a fermion self-energy at finite temperature. In Sec. Ill, the numerical results 
of the fermion spectral properties are shown. In Sec. IV, we discuss the gauge dependence of fermion pole appearing 
when T ^ m in an analytic way. Section V is devoted to a summary and concluding remarks. In Appendix A, we 
briefly describe how the abelian Higgs model is reduced to the massive gauge theory in the Stueckelberg formalism. 
In Appendix B, we present detailed calculational procedures for the fermion self-energy in our model. Appendix C is 
devoted to making an order estimate of some terms appearing in the text. 



II. U(l) GAUGE THEORY WITH MASSIVE GAUGE BOSON 

In this section, we formulate the U{1) gauge theory with a massive gauge boson, and introduce a propagator and 
a spectral function at finite temperature in the imaginary time formalism [13, 14]. We perform a calculation of the 
self-energy of a fermion coupled with a massive vector boson at one-loop order. 



A. General formalism 



First, we introduce a U{1) gauge theory with a massive gauge boson. The gauge boson acquires a mass by the 
Higgs mechanism, keeping the gauge symmetry. The gauge theory is one way to construct a renormalizable quantum 
field theory with a massive vector boson. We employ the Stueckelberg formalism [10, 11] proposed long ago, which is 
equivalent to the abelian Higgs model with a constant absolute value of the Higgs field [11, 15]. This correspondence 
is reviewed in Appendix A. Then our Lagrangian reads 



C = --F^.F^'' + (4^ _ ^ + _ igA^)^)^ + Cgf, (2.1) 



where A^, B and ip are a massive vector, a scalar and a fermion field, respectively. The scalar field B is called the 
Stueckelberg field, which corresponds to the phase of the Higgs field in the abelian Higgs model. F^y = d^A^ — d^A^ 
is a field strength, g the coupling constant, m the vector boson mass, and a is a gauge parameter. £gf is the gauge 
fixing term defined by 

Cgf = -^{^^.A^ + cxmBf. (2.2) 
la 

We work with the Minkowski metric, = diag (1, —1, —1, —1). We shall deal with a massless fermion assuming 
that the mass is neglected, which should be valid at high temperatures. Our Lagrangian is invariant under the gauge 
transformation except for the gauge fixing term, £gf: 

i,{x) e'^'^(^V(a;), (2.3) 
A^{x) ^ A^{x) + 9^A(a;), (2.4) 
B{x)^ B(x)+mk{x). (2.5) 

There are no interaction between the Stueckelberg field and the fermion field, and we chose the gauge fixing term so 
that the interaction term between the vector field and the Stueckelberg field vanishes. We can drop the Stueckelberg 
field as long as a correlation function is concerned, while it can not be when the thermodynamic potential is considered, 
where it is important to take into account the correct degrees of freedom. 



3 



p-k 




FIG. 1: The diagram which contributes to the fermion self-energy at one-loop order. The solid line represents the fermion and 
the wavy line represents the massive vector boson. 



The propagator of the free massive vector boson is now given by 

D,Ap) = (9,. - (1 - a)) . (2.6) 

In the a — > 00 limit, the propagator tends to 

D,Ap) ^ (9,. - ^) , (2.7) 

which is the massive vector-boson propagator in the Proca formalism^. 

The fermion propagator G(j)) in the imaginary time formalism [13, 14] is expressed with the self-energy S(p) as 

= j^r ^'-'^ 

where p° = iLJm = i{27Ji + 1)ttT is the Matsubara frequency for fermion. Note that G{p) and S(p) are 4x4 matrices 
with the spinor indices. The retarded fermion propagator is given by an analytic continuation, ia;,„ — > w + ie: 

G^ip,u:)^G{p,uj + ie)^— ^-—^ -, (2.9) 

^7 ^ P ■ 7 ^ ^ (Pi <-^) 

where the retarded self-energy is given by 

E«(p,cj) = E(p,cj + ie). (2.10) 

Introducing the projection operator on the (anti-)particlc sector A-|-(k) = (1 ± 7^7 • k)/2, we can decompose the 
retarded propagator and self-energy into the respective sector as follows: 

G«(p,c.) - G+(p,w)A+(p)7" + G-(p,L.)A_(p)7°, (2.11) 
Ei^(p, Lu) = E+(p, c^)A+(p)7" + S-(P, ^)A_ (p)7", (2.12) 

with E±(p,w) = Tt{T.^{p,uj)A±{p)j°)/2. 

In the particle sector, the pole Up = w(p) satisfies the following equation: 

G;i(p,a.p) = - IpI - E+(p,t^p) = 0. (2.13) 

From the analyticity of the retarded propagator, the pole is located on the real axis or the lower half-plane of complex 
w. If the imaginary part of the pole is small, the pole is well described in terms of a quasi-particle picture, where the 
real part of the pole corresponds to the energy while the imaginary part to the decay width of the quasi-particle. If 
the imaginary part is large, then it would be meaningless to consider excitations in terms of any particle picture. It 



^ Here we note that the propagator in the Proca formalism does not vanish but rather approaches a constant value in the p — ^ 00 
limit, in contrast to that in the Stueckelberg formaUsm. This causes the non-renormalizability and leads to a bad behavior at high 
temperature [8, 16, 17]. 
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is known that the self-energy at zero momentum has the following symmetry, 

ReE+(0,-w) = -ReS+(0,w), (2.14) 
ImE+(0, -w) =ImS+(0,cj), (2.15) 

which implies that if there exists a fermion pole at z = w — at zero momentum, there is also a pole at —uj — iT at 
zero momentum. 

Once the self-energy I]±(p,aj) is obtained, the spectral function of the (anti-) particle sector is expressed as 

1 T ^ / N 1 ImS±(p,a;) 

P± P,^ = — p,a.H— I ^2', T- 2.16 

TT TT (w =F IpI - R.ci;±(p, w))"^ + ImEj_(p,a;) 

When the peak is narrow enough, the position of the peak is given by Wp = |p| -I- Re S-|-(p, Wp) and the width of the 
peak is given by — Im I]+(p, Wp)/(2wp). 



B. Calculation at one-loop order 

Now let us evaluate the self-energy S(p) at one-loop order; the corresponding diagram is shown in Fig. 1. S(p) is 
expressed as 

(2.17) 



where Go is the propagator of the free fermion, I = p — k, and k'^ = iujn ~ i{2n + 1)T. Some manipulations lead to 

2 

E(p, iw,„) = -25^7^i?„(p, iw„i; m) H iB{p,iUjn; \/am) - B{p,iuJm]rn)] 

m^V \ \ J (2.18) 



(ai3(p, iw,„; ^/am) - i?(p, iw,„; m) ) ) -p'-f^ ( i?p(p, iw^; \/am) - Bf,{p,ium]m) 



The retarded self-energy I]^(p,a;) in the one- loop approximation is given by the analytic continuation iwm uj + ie 
from Y,{p,iuJm)- Here, we have introduced the following loop functions: 

B(p,.c.„.;^).tX: / (g3 (fc_^)2_^. ^- (2-19) 

B'^(p,.c.™;^).T^ j (^ (,„^)2„^2 fc2- (2-20) 

We see that there are two kinds of mass in Eq. (2.18), m and y^m, the latter of which is unphysical because it depends 
on the gauge parameter. However, the existence of such an unphysical mass causes two different high temperature 
limit as will be shown in Sec. IV. We will also show that S(p) approaches the fermion self-energy in QED if we take 
the massless limit to — )■ 0, which is not the case in the Proca formalism. It should be noted here that Eq. (2.18) shows 
that there is a special value of a: when a = 1, the terms containing the unphysical mass are all cancelled out and 
only the first term remains, i.e., Yi(ji,iuim) a=i = —2g''^^B^{p,iujm]'m). 

The self-energy in the (anti-)particle sector in the one-loop approximation is given by I]±(p,w) = 
Tr (S]'^ (p, a;)A±(p)7°)/2. Then, as is derived in Appendix B, we have for the imaginary part of S+(p,aj), 



ImS+(p,a.) = - f ^ I ^ dEf{f{Ef)^n{Ef-u:))[{~p'+m'a){\p\-u:Y + 2p'Ef{u:-\p\ 

e{~p^)[p\u:-\p\)TT^T'+L,{io-\p\){-ojp'-{oj-\p\){~p^ + m'a))] 



2 rE', 



f 



327r|p|2TO2 



' I dEfifiEf) + n{Ef - c^))[( V + m''){{\p\ - uof - 2w?) + 2[j? - 2m')Ef{io - |p| 



327r|p|2TO2 

,0 il—^^^P'^^^P' - 2m^)(^ - |p|)^^T2 +^.[2777^ -p^(|p|(a. - |p|) +to2)]], 
d27r|ppTO^ 

(2.21) 
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where 



IpI', Ef 



(^2 - |p|2 - m2)/(2(w ± IpD) and E'^ = (, 



am2)/(2(o.±|p|)). 



The real part ReI]+(p,Ct') may be obtained using the dispersion relation from the imaginary part. Especially, the 
finite temperature part of the real part of the self-energy, I]+(p,w)t/o = ^-(-(p,^) ^ '^+{v,^)t=o, is expressed as 



ReS+(p,a;)T^o 



1 

--P 

TT 



,ImE+(p, w')t^o 



(2.22) 



Here P denotes the principal value. The zero temperature part of Rel]+(p, lj) is not determined by Eq. (2.22) because 
it has ultraviolet divergence. We make renormalization using twice-subtracted dispersion relation, which reads 



ReE+(p,w)T=o = Co + ci(w - |p|) + 



(^-IpI)' 



dz 



ImS+(p, z)t=o 



We impose the on-shell renormalization condition, E+(p,a; = |p|) = and c?S+(p, a;)/9a;|^^|p 
and ci. The vacuum part ofImE+(p,aj) is obtained by taking the T limit of Eq. (2.21); 



(2.23) 

0, to determine cq 



Imi;+(p, w)t=o = 
Thus we arrive at 
ReI]+(p,a;)T=o 



g2 sgn(cj) 
327rm2 



(^-IpI) 



e{p^ - am^){p^ - am^Y - e{p 



2^{p^ + 2m2)(p2 _ m'^Y 



'■y- 



(2.24) 



? 9 

9 P 
'327r2m2 



(^-IpI) 



2m4 



(2 + a)- 



In 



-2a+ (3 + a2)- 



p2 - 



p2 — rn?a 



rn^a 



In 



a(p2 



p2 _ j7^2q. 



In 



(2.25) 



As mentioned before, our theory based on the Stueckelberg formalism approaches QED at high enough temperature 
where the masses are negligible in comparison with T. Let us see this. For T ^ gT ^ to, ^/am, the imaginary and 
real part of the self-energy are reduced to 



Imi;+(p,w)T-^o 



327r|p|2TO2 



0{-p')[p\^ - |p|)7r2r2 _ (p^ _ 2m2)(c. - |p|)^2j..] 



Rel]+(p,a;) 



T-i-oo 



16|p|2 

9^ 
16|p|- 



^r2(^-|p|). 



2|p| + (|p|-c^)ln 



(2.26) 



(2.27) 



respectively. Here, we have retained only the terms which are proportional to T"^ in Eq. (2.21). These Eqs. (2.26) 
and (2.27) coincide exactly with the well-known results in the HTL approximation in QED [1-4]. There is a caveat 
in the above manipulation, which has been taken for granted in the usual derivation of the HTL approximation in 
the gauge theory: The ignored terms may become comparable to terms which are proportional to in some gauges 
and hence the above naive power-counting turns out to be invalid. We will analyze this possibility in Sec. IV. 



III. NUMERICAL RESULTS 



In this section, we show numerical results of the fermion spectral function and the fermion poles at various tem- 
peratures. In the following, the coupling constant is fixed to a small value, g = 0.5, so that the analysis based on 
the one-loop calculation can be valid: Except when the coupling constant dependence of the pole is analyzed, the 
coupling constant will be fixed to g = 0.5. On the other hand, the gauge parameter, a, will be varied freely in order 
to see the gauge-dependence of the spectral properties of the fermion calculated at one-loop level. 



A. Low temperature (T <^ m) 



In this subsection, we show numerical results at a so low temperature that T dependence of the results is hardly 
seen, which may check our analytical and numerical calculations. 
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FIG. 2: The fermion spectral function in the particle sector p_|_ as a function of energy cu and momentum p for T = 0.4 m, g = 
0.5, a = l. 



Figure 2 shows the fermion spectral function in the particle sector (with a positive particle number) at T = 0.4 m 
for a = 1. There appears a very narrow peak near tu = |p|, which is very reminiscent of zero temperature case. This 
is natural for T <^ m, because the thermal effect is exponentially suppressed by the Boltzmann factor ~ exp{—m/T), 
and hence the breaking of Lorentz symmetry is small. This small breaking of Lorentz symmetry implies that the 
particle pole is almost on-shcll value at T = 0, i.e., lu = |p|, and hence the gauge dependence of the pole hardly 
appears. 



We plot the spectral function of the particle sector at T/m = 1.0, 1.5, 2.0, 2.5 for a = 1 in Fig. 3. We can see 
that the spectral function at these temperatures have structures qualitatively different from that at low temperature: 
Even at T/m ~ 1.0, we see a split of a peak around the origin seen in Fig. 2 into two peaks with a small bump in 
the negative energy region, which is reminiscent of the anti-plasmino peak known in QED/QCD at high T. These 
features are enhanced as T is raised, and we see a clear three-peak structure around |p| = at T/m = 2.0 with a 
prominent peak and a clear bump in the positive and negative energy region, respectively. We also see that the peak 
near the origin is attenuated as T is further raised up to 2.5 m. 

Since it is known that the details of the shape of the spectral function may be gauge-dependent in general, let us 
see how the three-peak structure depends on the gauge parameter. Figure 4 shows the gauge-parameter dependence 
of the fermion spectral function in the particle sector at T = 2.0 m; the gauge parameter is varied as a = 0.1, 1, 10 
and 100. One might find only single curve of the spectral function in the figure, although this figure actually shows 
four curves of it with different a; thus it clearly tells us that the shape of the spectral function at T = 2 m with a 
three-peak structure is virtually independent of the gauge parameter. 

The virtual gauge-independence of the shape of the spectral function implies that the pole of the propagator is 
also the case. We show the gauge (in)depcndence of the pole in the positive energy region at |p| = with a particle 
number in the left panel of Fig. 5, which shows that the pole position is almost independent of the choice of the gauge 
parameter, as anticipated: Note that the gauge parameter is varied in a wider range than in Fig. 4, i.e., a = 0.01, 0.1, 
1, 10, 100, 1000. A remark is in order here: The pole in the negative energy region at |p| = has the same properties 
as that in the positive energy region, as is assured by Eqs. (2.14) and (2.15). 

Such a gauge-independence of the poles necessarily reflects in that of the spectral function. The right panel of 
Fig. 5 shows the fermion spectral function at zero momentum for the wide range of a up to 1000, together with 
that obtained in the Proca formalism'^. From this figure, we confirm that the spectral function at zero momentum is 
virtually gauge-independent for the wide range of a. 



3 The rapid decrease of the spectral function in \ui\ < 0.1 is caused by the exponential damping of Im S-|-(0, cj) in that region. 



B. Intermediate temperature (T 



m) 
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FIG. 3: The fermion spectral function in the particle sector p+ as a function of energy uj and momentum p at <; = 0.5, a = 1 
at various temperatures. 




FIG. 4: The spectral fermion function in the particle sector p+ as a function of energy u and momentum p at T = 2.0 m, g = 0.5. 



We also note that the position and the width of the peaks coincide with the real and imaginary part of the poles, 
respectively, which is due to the fact that the imaginary part of the poles is small in comparison with the real part, as 
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FIG. 5: The left panel shows the gauge dependence of the fermion pole with a positive particle number in the positive energy 
region for T = 2.0 m, g = 0.5 at zero momentum: the horizontal and vertical axes denote the real and imaginary part of the 
energy, respectively. The right panel shows the fermion spectral function in the particle sector p+ as a function of energy oj for 
the same T, g as the left panel. 




FIG. 6: The left panel shows the coupling constant dependence of the fermion pole with a particle number in the positive energy 
region for T/m = 2.0 and a = 1.0, at zero momentum: the horizontal and the vertical axis denote the real and imaginary 
part of the energy, respectively. The dots from left to right correspond to g = 0.1, 0.2, 0.3, 0.4, 0.5, respectively. The right 
panel shows the coupling constant dependence of the real part of the fermion (particle) pole in the positive energy region at 
zero momentum for the same T/m and oi. The solid line is the fitted linear function. 



seen in the left panel of Fig. 5. Thus the shape of the spectral function with a three-peak structure necessarily gets 
to have almost no gauge-dependence. 

We show the coupling constant dependence of the fermion pole at zero momentum in Fig. 6 for T/m = 2.0. The 
real part is almost proportional to g, like that in QED in the HTL approximation. The coupling constant dependence 
of the imaginary part is not large. 

What is the mechanism for realizing the three-peak structure of the fermion spectral function? Figure 7 shows the 
real and imaginary part of the self-energy for T/m = 2.0 and a = 1.0 at |p| = 0, together with the corresponding 
spectral function. A detailed analysis of the imaginary part tells us that the peaks of the imaginary part correspond 
to a Landau damping of the fermion by a scattering with thermally excited bosons. Since these features of the fermion 
self-energy is very similar to that shown in [8], the mechanism for realizing the three-peak structure found in our 
formalism is understood to be the same as discussed in [8]. 
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FIG. 7: The upper figure shows the fermion spectral function in the particle sector p+ as a function of energy u for g = 0.5, T = 
2.0m, and a = 1. The lower figure shows the real part and the imaginary part of the fermion self-energy as functions of energy 
bj for the same g, T, and a. The dotted line denotes oj. The intersection points between this line and ReE+(0,a;) are the 
solutions of the dispersion relation. 



C. High temperature (T ^ m) 

In this subsection, we show numerical results in the high temperature (T ^ m) region, where the mass of the vector 
boson (and the fermion) can be neglected in comparison with T, i.e., m/T — > 0; this means that T itself may not be 
infinitely large. 

We show the fermion pole in the positive energy region in the left panel of Fig. 8 at T = 40.0 m and for a = 0.1, 1, 10 
and 100. The pole in the Proca formalism and that in the HTL approximation in QED are also shown. We see that 
the gauge dependence of the fermion pole is no longer negligible. Since the exact pole position in the complex energy 
plane should be gauge-independent [12], the above result suggests that the one-loop analysis is no longer valid in 
this high-T region in contrast to the T <^ m and T ^ m regions, at least in some gauge. We will present a detailed 
discussion on how the gauge dependence arises at high T region in Sec. IV. 

One should notice that the pole for a = 0.1 is located in the upper energy plane, which could be problematic 
because it implies a loss of the analyticity of the retarded propagator and also negativeness of the spectral function, 
as seen from Eq. (2.16). 

For a ~ 1, there appear clear two peaks in the spectral function in a robust way, as shown in Fig. 9; the two peaks 
are found to tend to the normal fermion(particle) and the anti-plasmino of QED in the HTL approximation [1-4], 
respectively; see Sec. IV. 

There persists the other peak at the origin in the energy-momentum space. One can confirm that its residue is 
of the order of m'^ / {g'^T'^) , which is very small if we consider the gT ^ m case, by making power counting. Such 
a peak at the origin was also obtained in [8], though in the Proca formalism. One should also remark that such a 
peak at the zero energy is not obtained in QED in the HTL approximation, in which the vector boson mass is set 
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FIG. 8: The left panel shows the gauge dependence of the fermion pole at p = in the positive energy region at T = 40.0m for 
g = 0.5 : The horizontal and the vertical axis denote the real and imaginary part of energy. The right panel shows the fermion 
spectral function in the particle sector p+ as a function of energy uj for the same T and g. 
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FIG. 9: The fermion spectral function in the particle sector p+ as a function of energy cj and momentum p at T = 40.0m with 
g = 0.5 and a = 1. 



to zero from the beginning, in contrast to the present case^. It should be intriguing to explore whether this peak at 
the origin extends to a finite- |p| region, and hence the three-peak structure of the fermion spectral function persists 
even in such a high-T region, i.e., for T 3> m. In fact, this is a challenging problem in quantum field theory at finite 
temperature, because a sensible analysis of such an infrared region requires a systematic method to remove the so 
called pinch singularities [18]. This task is beyond the scope of the present work, and we leave such an analysis as a 
future work [19]. 

Our numerical calculation has shown that one can have virtually gauge-independent results even in the one-loop 
analysis if the gauge parameter is in the region a 1. We shall argue that the pcrturbative expansion should be 
valid for a <^ 1/g in Sec. IV. It means that the spectral function of a fermion coupled with a massive vector boson 
as calculated in the Stueckelbcrg formalism nicely approaches that in QED in the HTL approximation at high T 



The absence of a peak at the vanishing energy in QED with the HTL approximation is easily understood as follows: In the HTL 
approximation of QED, ReS-j_(0,[^) behaves as l/uj in the ui limit. Thus at ui = 0, the pole condition Eq. (2.13) will not be satisfied 
and hence there can not exist a pole at the origin. 
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irrespective of the choice of the gauge parameter a, if the order of a is confined to a ^ 1/ g. This is actuahy aheady 
suggested by the asymptotic form Eq. (2.27) for T ^ m, ^fam. 

IV. ANALYSIS OF GAUGE DEPENDENCE OF THE POLE AT HIGH TEMPERATURE (T > m) 

Our numerical calculation has shown that the pole position of the fermion propagator is virtually independent of the 
gauge parameter a for the cases of T <C m and T ^ m: The former case is simply because the thermal contribution 
due to a boson with a mass m is greatly suppressed by a Boltzmann factor exp(— jti/T) <C 1 when T <^ m. By 
contrast, for T ^ m, the numerical results in Sec. IIIC show that the pole of the fermion propagator has a large 
gauge dependence for large a. In this section, we discuss the gauge dependence of the pole of the fermion propagator 
at weak coupling at high temperature. In particular, we focus on the region gT ^ m. In this region, one expect that 
the mass of the vector boson can be neglected, and thus the self-energy approach that in the HTL approximation of 
QED [1-4], in which the fermion has the pole of order gT. Therefore, we analyze the pole of the fermion propagator 
by assuming lu ~ gT. Here we introduce a small dimensionless parameter, 

A^(^)«l. (4.1) 

Thus we have two small dimensionless parameters, g and A, which are treated as independent parameters, so that 
the self-energy is expanded by combined powers of g and A. If the the power of g and A are both positive, the high 
temperature limit will be well defined and smoothly connected to that of QED. However, as will be shown below, an 
inverse power of A appears at one loop level when the gauge parameter is large, and hence the high temperature limit 
becomes inevitably different from that of QED. 

In the following analysis, we put |p| = for simplicity. The pole position obtained in the perturbation theory 
generally depends on the gauge parameter as well as g, A and T due to the truncation of the perturbative expansion. 
We parametrize the pole Wpoie of the fermion propagator as 

Wpoic-5rF(5,A,a), (4.2) 

where F[g, A, a) is a function of order one, and depends on the gauge parameter a. If the limit, 

Fo= limf^(5,A,a), (4.3) 

is independent of a, then the pole is independent of the gauge parameter at the order gT. Thus one sees that the 
gauge dependent part may be defined by 

(JwpoioCg, A, a) = ujpoicig, A, a) - Wpoio(A), (4.4) 

where wOqi^IA) = gTFo{X). For a reference, we recaU that Fo(A = 0) = 1/(2V^) in the case of QED [4]. When the 
inequality, 

'^poic(A) > Sojpo\c{g, A, a), (4.5) 

is satisfied, the gauge dependence can be neglected. In reality with a finite g, the region of the gauge parameter 
satisfying Eq. (4.5) will be limited. We shall call the region that the gauge parameter satisfies Eq. (4.5) as an 
adequate gauge parameter region. The purpose of this section is to find the adequate gauge parameter region. 

Let us first show a numerical result of the real and the imaginary part of the pole at T = 40.0 to as functions of 
a in Fig. 10: For a large a (1 <C a), the a dependence of the real part of the pole is large, and especially for very 
large a, say a '--^ 3 x 10'*, the magnitude of it is no longer of 0{gT), but is of a smaller order, 0{m), as will be shown 
later. The imaginary part of the pole for a ^ 1 is positive and apparently problematic because it means that the 
analyticity of the retarded propagator is lost and the fermion spectral function will become negative. As we shall 
show later, however, the absolute value of the imaginary part is of 0{g^T) and should be considered together with 
higher order contributions. So the negative imaginary part with a small absolute value can be ignored in this order 
of the coupling. 

Now we shall show that such an order estimate of the pole can be done analytically. We start with an analysis of 
the self-energy, under the condition that ui ~ gT, by decomposing the self-energy (2.18) to seven parts, 
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n=l 



(4.6) 
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where we have introduced the following dimensionless functions: 





2 

= —B°{0,Lo;m), 


(4.7) 




2 2 

= H ^i>(0, y'am), 


(4.8) 




= 2 5(0, w;m), 


(4.9) 


C4(w) 




(4.10) 




= +.9^5(0, a;;m), 


(4.11) 






(4.12) 




= +^B°(0,tj;m). 


(4.13) 



Here B(0, w; m) and B"{0 ,Lu; m) are obtained by performing the analytic continuation {iujm ^ uj + ie) to B{0, iuJm] fn) 
and i?°(0, zw™; m). From Eq. (2.13), one sees that the poles satisfy the condition 

C(Wpole) = 1. (4.14) 



We first note that if the following equation has a root, 

limC(5TF) = l, (4.15) 

there is a pole Wpoie = gTF of order gT. Furthermore, if the F happens to be independent of a, the pole is gauge 
independent. Therefore, let us take the following function as a measure of the gauge dependence, instead of Eq. (4.4): 

5C(g,A,a)sC«J-l. (4.16) 

Then a criterion for the adequate gauge parameter region may be given by 

<5C(5,A,a) < 1. (4.17) 

We now make an order estimate of the seven terms defined in Eqs. (4.7) ^ (4.13). One finds that this task is 
reduced to that of B{0, uj; m) and i?°(0, a;; to). The following relations are shown in Appendix C: 



When T > w > M, 



When T > M > w. 



B{0,u:-M)^-, ^ ' ' (4-18) 



, r ^ for Af 2 < ujT 
B(0,u}:M) ^ < %.2 

^ ' \^ for A/2 > wT 

SO(0,w;Ar) / 4 for AP ^ uT 



(4.19) 



for A/2 > LoT 



• When M:^T:^uj, 

,<„,.M,.I1, ^^^^ 

Here, A/ denotes m or ^/am, and the vacuum parts have been dropped. Using these relations, we will find the 
adequate gauge parameter region in the following subsections, and the result is summarized in Fig. 11, which shows 
that the adequate gauge parameter is a <^ I/5. 
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20 30000 
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FIG. 10: The upper panel shows the gauge dependence of the real part of the pole in the positive energy region at zero 
momentum at T = 40.0m, for g = 0.5. The horizontal and the vertical axes denote the gauge parameter a and the real part of 
the pole, respectively. The data whose a is smaller than the upper limit of adequate gauge parameter, 1/g = 2.0, are denoted 
by the solid line and the data whose a is larger than 1/g by the dotted line. The data which satisfy a > l/(g^A) — 1600, which 
are expected to approach the value in the Proca formalism, are plotted with another dotted line. The upper and lower dashed 
lines show the thermal mass in the HTL approximation gT/ {2^/2) and that in the Proca formalism \/6m/Vl + 48 A ~ 2.31m, 
respectively. The lower panel shows the gauge dependence of the imaginary part of the pole for the same T and g as in the 
upper panel. 



Ci(a;), (73(0;), C5(w), and C-;{uj), which do not depend on the gauge parameter, are estimated to be 

CaM^ffQ), CsH^g, C,{lo) ^ . (4.21) 

We remark that Ci{lo) coincides with the fermion self-energy in the HTL approximation in QED, as it should be. 

We note that 6*3(0;) and 6*7(0;) are of the order of an inverse power of A, which would make it impossible to take 
the massless limit. These 'dangerous' terms are found to be nicely canceled out with other Ci(o;)'s when a is not so 
large, whereas for large a, the cancellation does not happen, and the condition 5C ^ 1 can not be satisfied, as will 
be shown below. 



A. 1 2> aA (case 1) 

When 1 ^ aA, the mass scale m and y/am are negligible in comparison with (7T, and then the self-energy coincides 
with that in HTL approximation in QED in this approximation, because the non-leading terms can be neglected as 
we have seen in Sec. II. 

On account of the order estimate Eq. (4.18), we obtain 

C2(w) ~ 9 , Ci{Lu) ^ ga, Ce{oj) ^ . (4.22) 
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(a) g>X 



1 







1/g 1/X l/(g^X) oo 

a 



(b) g<X 







l/X l/g 



a 



l/{g'X) oo 



FIG. 11: The classification of tfie gauge parameter regions depending on the relative magnitude of A = w? /{gT)"^ to g at high 
temperature, i.e., for gT ^ m or A ^ 1. Each number in the figures corresponds to the case with same number in the text. 
The light gray regions are the adequate gauge parameter regions for which Eq. (4.5) is satisfied, whereas the dark gray region 
are not. 



The leading terms of C2{(jj) and Cq{ijj) cancel out^ with Cz{oj) and Ct{lS), respectively, and the terms of the order 
ga and g remain^ . Thus one sees that the gauge dependent part is of the order of ga, which means that the pole 
is gauge independent in practice, which also can be confirmed from Fig. 10, provided that the inequality a ^ l/g is 
satisfied; in this case, the adequate gauge parameter region is the region of a which satisfies the above inequality. 

Let us see that the imaginary part of the self-energy is of order g^T in the Landau gauge (a = 0). For a = 0, the 
imaginary part of the self-energy is evaluated to be 



lmi;+(0,a;) 
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2 r / , .2 , .2 



647rm2 

lOTT 



6.3 



tanh h coth - -^^ — coth h tanh 



(4.23) 



where the first line is obtained by substituting p = and a = in Eq. (B2). In the second line the inequality 
T 3> w 3> m was used. Although this is positive and apparently breaks the analyticity, this order of the coupling 
should not be determined in the one-loop order, because the two-loop diagrams contain contributions of order g^T . 
In fact, the analyticity problem can be cured by taking into account the two-loop diagrams. 

B. 1/p^ > qA > 1 (case 2) 

1. Q < l/(ffA) (case 2-a) 
On account of the order estimate for C2(w), 6*4(0;), and Cq{uj) given in Eq. (4.19), we obtain 

C2{u) ^ g , Ci{u) ^ ga, C^{u) ^ , (4.24) 



^ As seen from Eqs. (C14) and (C19), B{p;m) and B^{p;m) yield terms which are proportional to m, so C2(aj), C3{u}), Cs(lu), and C7(aj) 
seems to yield terms which are proportional to and the massless limit {m/T — >■ 0) can not be taken. Actually, from Eqs. (C14) 

and (C19), we see that the terms discussed above cancel out. Therefore we can take the massless limit and the fcrmion self-energy 
approaches that in QED. 

® These terms come from the imaginary part of B{0, uj; m), -8*^(0, uj; m), -6(0, uj; y/am), and -B"(0, u): y/am), which can be confirmed by 
retaining the next-to-leading term in Eqs. (C7) and (C15). The contribution from the real parts arc much smaller than these terms. 
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which is the same as Eq. (4.22). The leading term of Ce{uj) cancels out with that of Cr{uj) and the remaining terms 
are of the order of ga^ , while C2(w) and 6*3(0;) do not. In the present case, however, 64(0;) is larger than 6*2(0;) 
and 63(0;), and hence 6*4(0;) dominates the gauge-dependent terms, which should be made small. This smallncss is 
guaranteed when a <C 1/g, which defines the adequate gauge parameter region; this coincides with that in the case 1. 
Notice that the inequality g < X must be assumed in this case, otherwise the inequality a < 1/g can not be satisfied. 



2. Q > l/igX) (case 2-h) 

Again, on account of the order estimate in Eq. (4.19), 62(0;), 64(0;), and 65(0;) are estimated to be 

Though 64(0;) and 67(0;) have the same order of magnitude, they do not cancel out, which can be confirmed from 
Eqs. (C20) and (C26). Therefore the largest contribution is 64(0;) ^ 1/A 1; the present region for the gauge 
parameter (a 3> l/(.gA)) is not an adequate gauge parameter region. 



C. a\ » l/g^ (case 3) 

Here we treat the case where the gauge parameter is far larger than 1/Xg^; this case includes the unitary gauge 
{a — >■ 00). Owing to the order estimate Eq. (4.20) for this case, 62(0;), 64(0;), and 66(0;) are estimated to be 

1 /1\' „ , . /1\ „ . . 1 



62(.)^-^-J , 64(.)^^-j, C^i^)-^. (4.26) 

Again, although 64(0;) and 67(0;) have the same order of magnitude, they do not actually cancel out on account of 
the difference in the coefficients, which can be seen from Eqs. (C20) and (C42). In this case, the largest contribution 
is 64(0;) ^ 1/A ^ 1; therefore this region is not an adequate gauge parameter region. This also suggests that the 
order of the pole in the unitary gauge is not gT. 

So let us now discuss the pole in the unitary gauge. First we assume the pole is of order m instead of gT. Using 
Eqs. (4.20), (C50), and (C54), we have the following order estimates: 

6.(.)^i, C2(.)^-^, ^^3(0;)^-^, 64(0;)^, C,i.)^j=, 

(4.27) 

For a — > 00, 62(0;) and 65(0;) vanish. Furthermore 63(0;) and 65(0;) can be also neglected, because gVX = m/T <C 1. 
The remaining parts 61(0;), 64(0;) and 67(0;) are estimated more precisely with the use of Eqs. (C20) and (C42), as 
follows: 

^'M=^8xS> '^'(->==4 ^'M=^-T^' 



Collecting these terms, we reach at 



This is precisely the same as the result obtained in the Proca formalism in the high temperature limit [8] ; the thermal 
mass in this case reads 



VTTi8X 



V6m, (4.30) 



which can be seen from Fig. 10 also. We have confirmed numerically that the spectral function also approaches that 
in the Proca formalism as a — >■ c». 



This term comes from lmB'^{0,uj; y/am), which can be confirmed from Eq. (C30). The real parts are negligible compared with this 
term. 
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D. Brief Summary 

Let us summarize and diseuss the results obtained so far in the preceding subsections for the gauge-parameter 
dependence of the fermion propagator at high temperature. The results are summarized in Fig. 11. We have found an 
adequate gauge parameter region as a ^ l/g in which possible gauge dependence is of higher order of the couplings 
and hence can be neglected. We remark that this parameter restriction should also apply to QED. 

In fact, the electron self-energy in the HTL approximation in QED at next-to-leading level in the one- loop order 
reads [20, 21] 

^"c - + &y (-'4 - + <' - (-1° ^ ^) ■ 

If it were that uj ^ eT and a ~ 1/e, the third term would be the same order as the leading term, which implies that the 
gauge-independence is badly broken. Conversely speaking, if a is much smaller than 1/e, then the gauge-dependent 
part becomes of order and can be neglected; this gives an adequate gauge-parameter region for QED. 

There is a difference between QED and the massive vector theory in the Stueckelberg formalism: The limit a — > oo 
can not be taken in the former case because the third term in Eq. (4.31) diverges, while it can in the latter case. In 
this limit, the pole of the fermion propagator in the latter case becomes of order m <^ gT . 

On the other hand, also small a yields the problem; the imaginary part of the fermion self-energy becomes positive 
for a ^ 1, which breaks the analyticity of the self-energy. 

The a dependence of the self-energy can be understood intuitively as follows: In the Stueckelberg formalism, there 
are two masses: the physical mass m and the unphysical mass ^/am. For gT ^ m, y/am, the self-energy naturally 
approaches that of QED because both masses can be neglected. By contrast, if the unphysical mass ^/am is not 
smaller than gT, ^foan can not be neglected, although the temperature is enough high compared with the physical 
mass TO. In this case, the self-energy at one loop level does not approach that of QED. 

We also note here that in the g <^ X case, the boson mass to can not be taken to zero from the outset, but the 
self-energy is approximately equal to that in the HTL approximation in QED as long as a <i; l/.g. This implies that 
the value, Tg"^/^, is the upper limit of the boson mass that we can neglect. 

Summarizing the situation, we see that there are three region of the gauge parameter: In the first region, the theory 
approaches QED, and is in adequate gauge parameter region. In the second region, the theory approaches QED, but 
is out of adequate gauge parameter region. In the third region, the theory does not approach QED, and is out of 
adequate gauge parameter region. We concludes that the gauge parameter should be chosen to be 1 < a ^ l/.g in 
numerical calculations in this formalism. 



V. SUMMARY AND CONCLUDING REMARKS 

We have investigated the spectral properties of a fermion coupled with a massive vector boson in the whole tem- 
perature (T) region at one-loop order. The vector boson with a mass (to) is introduced as a C/(l) gauge boson in the 
Stueckelberg formalism so that the high T limit, or equivalently the massless limit in the sense that m/T — )■ 0, can 
be taken^: We have successfully analyzed and clarified the characteristics of the spectral properties of the fermion in 
the distinct three regions of T, i.e. (I) T <C m, (II) T ^ m and (III) T ^ to regions, in a unified way. We have also 
carefully examined the possible gauge dependence of the spectral properties of the fermion in the respective three T 
regions, separately. 

In the region (I), the fermion spectral properties hardly change from those in the vacuum, which are gauge- 
independent. In the region (II), the fermion spectral function gets to have a three-peak structure in the small 
momentum region with supports in the positive, zero and negative energy regions; the three-peak structure becomes 
prominent when T ~ 2m for g — 0.5. We have confirmed numerically that the fermion poles and hence the fermion 
spectral function shows virtually no dependence on the gauge parameter (a) for T m. It is thus natural that the 
similar three-peak structure of the fermion spectral function was obtained in the Proca formalism for the massive 
vector field for T ~ m [8], since the Proca formalism exactly corresponds to the unitary gauge a — > cx3. Conversely 
speaking, the three-peak structure found in [8] is not an artifact by a special choice of the gauge and is physical. 

It is interesting that the spectral function of a fermion coupled with a scalar massive boson shows also a similar 
three-peak structure for T ^ m at the one-loop level [8, 9]. The present analysis has established that a fermion 



As is mentioned in Introduction, the Proca formalism does not yield sensible results for T 2> m [8, 16, 17]. 
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coupled with a massive boson with a mass m has a three-peak structure for small momenta with supports in the 
positive, vanishing and negative energy regions at temperatures comparable with the boson mass, irrespective of the 
type of the boson, at the one-loop order. 

For T ^ m (region (III)), the fermion spectral function tends to have distinct two peaks precisely corresponding to 
those seen in QED in the HTL approximation [1-4]. It means that our formalism nicely describes the spectral function 
of the fermion coupled with a massive vector boson even in the high T region. There is, however, a tricky point related 
to a possible gauge dependence. We have found that there exists an adequate region of the gauge parameter a in 
the high-T region for the perturbation theory at finite T: If a is of the order 1, the analysis at the one- loop order 
makes no problem and is reliable, keeping the positivity of the spectral function and so on; otherwise, however, these 
fundamental properties may be lost. This is because there exist two mass parameters, i.e., the vector boson mass m 
and the ghost mass ^/am inherently in the Stueckelberg formalism. Thus the precise high-T region should be defined 
by the two conditions, T ^ m and T 3> \/am. Our extensive analytic study has proved this observation and showed 
that when a l/g {g is the coupling constant), the one- loop analysis is reliable even in the region (III). Accordingly, 
if the unitary gauge (a — >■ oo) is adopted for the massive vector boson, the one-loop analysis can not be valid for 
T ^ m, as is shown in a different context [16, 17]. 

Our numerical calculation has shown that there still remains a peak at the origin in the a;-|p| plane, though with 
a faint strength even in the region (III); this is in contrast to QED in the HTL approximation where such a peak is 
absent. Although it is an interesting possibility that the three peak structure persists at the high-T region and even in 
QED, a sensible analysis of the spectral properties around such a low-energy region requires a resummed perturbation 
theory to deal with possible pinch singularities [18]. Thus we leave an analysis of the spectral properties in the very 
low energy region as a future work and hope to report elsewhere [19]. 

We can think of some physical situations where the present analysis can be relevant, since massive vector bosons 
at finite T appear in various physical systems. In QCD, vector bosons or vector-bosonic modes may decrease their 
masses in association with the restoration of chiral symmetry at finite T [22] . It would be not surprising if there exist 
vector bosonic modes even in the deconfined and chiral symmetric phase in the vicinity of the critical temperature 
Tc, since the existence of other hadronic modes [7, 23, 24] and bound states [25] are suggested in that temperature 
region. There may also exist a vector-type glue ball in such a system. Then the present analysis would suggest that 
the quark spectra can be largely affected in such a system where the boson mass is comparable to T in the order of 
magnitude. In the electro-weak theory, the dispersion relation of neutrinos at high T may possibly be affected by 
the weak bosons the masses of which change with T [26]. One of the findings of the present analysis tells us that 
the one- loop analysis in the unitary gauge can not be applicable when m/T <^ 1. which includes the vicinity of the 
critical point. 
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Appendix A: THE CORRESPONDENCE BETWEEN THE STUECKELBERG FORMALISM AND THE 

ABELIAN HIGGS MODEL 

In this Appendix, wc briefiy show that the abelian Higgs model is reduced to the U{1) gauge theory with a massive 
gauge boson in the Stueckelberg formalism [11. 15]. 
The Lagrangian of the abelian Higgs model reads 

Cmsss = -\f^^-F^, + \id^ - zeA^)*!^ (Al) 

with F^i, = dfjAi, — d^A^. Here and $ denote the vector and the Higgs field, respectively. 
We fix the absolute value of the Higgs field $ and use the following polar representation: 

m / ieB(x) \ , , ^ 

$ = — =exp ^ , A2 

eV2 \ m J 
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with |$| = m/(e-\/2). We remark that the scalar field B, which will turn to be identified with the Stueckelberg field, 
is introduced as the phase of $. Then, the Lagrangian becomes 

1 / 1 \^ 

/iHiggs = -^F^^F^. + — (^A^ - -a^Sj , (A3) 

which exactly gives the free Lagrangian of the massive vector field in the Stueckelberg formalism and its interaction 
term with the Stueckelberg field B given in Eq. (2.1) in the text. This is what we wanted to show. 



Appendix B: CALCULATION OF THE FERMION SELF-ENERGY 

Here we derive Eq. (2.21) in the text. We first recall that the retarded self-energy S]^(p,w) in the one-loop 
approximation is given from Eq. (2.18) by the analytic continuation, iaj,„ oj + ie = p"; see Eq.(2.10). Then its 
imaginary part reads 



m 



2 

TO^ L \ \ 

2 



Iml]'^(p, w) = - 2g^jnmB^{p,w;m) + ^ \ ^(p^ (lmB(p,w; V^m) - Imi?(p 



(Bl) 

(aIm_B(p, uj; y/am) — Imi3(p,a;; m) J J — p^7^ (lm_B^,(p, lo; ^fam) — Imi3p(p, m) 



where = cj^ — |pp. By applying the projection operator for the particle sector, we have 
ImS+(p,L.) =iTr[ImI]«(p,c^)A+(p)7°] 



2 r- 

9 ' 



2.g2 (lmB°(p,a;;TO) - pJmBi(p, w; m)^ 

{lo — IpI) (^p^ ^Im_B(p, Lo; \farri) — Imil(p, a;; m)^ — ^almi3(p, a;; \foLm) — Im_B(p, w; rn)^ 
— p? ^Inii3°(p,a;; \fomi) — Im_B"(p,a;; jti) — pjmi3j;(p, cj; \/am) + pjmi3j(p,a;; m)j . 

(B2) 

Here wc have used the relation Tr[^A+(p)7'']/2 = — |p| and Tr[i3''(p, a;; TO)7pA-)-(p)7°] = i?"(p,aj;m) — 
PiBi{p,uj;m). Notice that lmB{p,uj;m) and ImB'^{p,uj;m) contain 6{sEf — tEi, — uj) with s,t = ±1, Ej = |k|, 
and Eb = ^ ni? + (p — k)^, while Imi3(p, w; -y/am) and Imi3''(p, cj; ^am) contain S{sEf — tEj^ — u) with E^ = 
■\/ am'^ + (p — k)2. Collecting the terms which contain 5{sEf — tEb — u), we have 

d^k v-^ stir 



{SisEf - tEb -u) term) ^g' J J2 j^J{sEf - tEi, - uj){f{sEf) + n{sEf ~ u)) 



X -2isEf - p • k) - —{u - \p\)ip' - m') + —p\sEf - p • k) 



(B3) 



Here we have used the following formulae; 

ImB(p,c.;m)=y'|^ ^§^J{sEf - tE^ - uj){f{sEf) + n{sEf - oo)), (B4) 

lmB^{p,uj;m)^J J2 -^^^5{sEj - tE^ - u:)if{sEj) + nisEf ~ u;))sEf , (B5) 

E ^^/S{sEf-tEb~u;){f{sEf)+n{sEf-u)). (B6) 
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Using the constraint p • k = i?/ cos9 = (m? — + 2LusEf)/2\p\ posed by the delta function, we have 

E -^^HsEf-tEt-u;){fisEf)+nisEf~u;)) 



s,t=±l 

-2-^.2 2m2-p2 



-iiu - IpI)^ - 2m') + sEf—^iu |p| 



2|p|m2 '^'^ ' ' to2|p| 



s,t=±l 

X (/(si;/) + n{sEf - w)) / dEbS{sEf - tE^b - w) 



(B7) 



with e± = {\p\ ± -E'/)2 + m2. Here the momentum integral is converted to that for the particle energies through 
the formula 

27r f°° 

(fk= -- / dSf^;/ / dEbEb. (B8) 
IpI Jo Je. 

We now evaluate the integral separately for the time-like and the space-like region. For the time-like region 
> |p|), we have 

{S{sEf - tEb - io) term) =^^^^ dEf{f{Ej) + n{Ej - c.))[( V + m'){{\p\ - - 2m') ^^^^ 

+ 2{p' -2m')Ef{uj-\p\)], 
with E'^ = (p2 — m2)/(2(cLi ± |p|)). For the space like region < |p|), we have 

92 r^7 

{S{sEf - tEb - u:) term) =^^-^__ J dEf{f{Ef) + n{Ef - uj)) 

X [(-p2 + m'){{\p\ - uj)' - 2m') + 2{p' - 2m')Ef{uj - |p|)] (BIO) 

~ 327r|p|2m2 [^^' " ^"^'^^^ " 'P')^'^' + " ^'^'P'^^ " 'P'^ + 

Here we have used the integral formulae 

rfi?/(/(i?/)+7l(i?/-C.))=-^, (Bll) 

DO 

dEfifiEf) + n{Ef - u))Ef =^ - ^. (B12) 



Combining Eqs. (B9) and (BIO), we have 

q' f^l 

{S{sEf - tEb - CO) term) ^ ^,J,, , / , dEf{f{Ef) + n{Ef - u)) 



327r|p|2m2 J 

X [(-p2 + m2)((|p| - a;)2 - 2m2) + 2{p' - 2m')Ef{uj - |p|)] (^13) 
e{-p')[{p' - 2m'){u - IpDtt'T' + Lo[2m' - p'{\p\{u - |p|) + m')]]. 



327r|p|2m2 
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Next, we collect the terms which contain 5{sEf — tEj^ — oj); 

q2 poo 

i6{sEf - tE', - w) term) ^ \ / dEj ^ st{u: - |p|)[-(c^ - |p|)(p2 _ am^) + 2p^sEf)] 

6ZTrm ipi ^^^^^^ 

X ifisEf) + n{sEf -uj)) [ ^ dEl5{sEf - tE^, - uj) 

J e'_ 

dEfifiEf) + n{Ef - c^))[( V + rn^a){\p\ - ujf + 2p^Ef{uj - |p 



32^|p|2m2 7^.+ 

-P')[p'{^ - \P\)7^'T' + ^{^ - |P|)(-V - - |P|)( V + rn'c.))], 



327r|p|2TO2 

(B14) 



where e'^ = ^ (|p| ± -B/)2 + Q:m2 and i?^ = (p2 — aTO2)/(2(cj± |p|)). Combining Eqs. (B13) and (B14), we arrive at 
Eq. (2.21). 

Appendix C: POWER COUNTING OF THE FERMION SELF-ENERGY 

In this appendix, we make an order estimate of B{p,uj;m) and B'^{p,uj;m) used in the analysis done in Sec. IV. 
The estimate will be made separately for the following four limiting cases: T ^ ui ^ T ^ m ^ uj , m ^ T ^ u 
and T ^ to ^ m. For simplicity, we set |p| =0. 

We start with _B(0,icjm;m) and B'^[0,iLOm]'rn) in the imaginary time formalism: 

i?~o(o,..,„;H ^[^T sE,nsE,) + itE, + ^.MtE,) 

Performing the analytic continuation and dropping T = part, we get the corresponding retarded functions, the real 
parts of which read 



^ E ^ AE,EfsEf-tE,-J ^^'''^^ + - ™)^=° 



(C3) 



27r 



/■ 5i sEff(sEf) + (tEb+uj)n(tEb) -r,, 



cj f°° f ^ 1 E'i 2E^f + m^ - u? 
:-^P / 2Ep{Ef)—^ ^ + n{Ei,)-f ^ ^ — 5- 



(C4) 



Here we have defined T ^ parts as 5(0, w; m)T^o = B{0,uj;m)~ B{0,uj; m)T=o and i?°(0, w; m)T^o = i3°(0,a;; m) — 
_B"(0, w; m)T=o- In the evaluation of these integrals, we have to carefully separate and deal with the contributions 
from the soft region [Ef <^T) and the hard region [Ef ^ T) unlike the calculation of the leading order of the HTL 
approximation, because the contribution from the soft region can be much larger than that from the hard region in 
some limits. 

From Eqs. (B4) and (B5), the finite temperature parts of Imi3(0, w; m) and lmi?°(0, w; m) are calculated to be 
ImB(0, c.; m)T^o =^ / dEf-f [~{f{Ef) + n{Eb))S{Ef - Eb + uj) - {f{Ef) - n{Eb))5{Ef + Eb - uj)] 

,2 / f w? -u'^ \ /™2 , . .2 \ \ y^'') 

sgn(m^ - w^)/ I I + n 



167ra;2 



2w 



m'^ + uj^ 



2u} 



21 



and 

1 E"^ 



1 r°° E 

" (C6) 



327ra;3 



2w 



2w 



respectively. In the first lines, the positiveness of a; has been taken into account. 

1. T > > m 

a. B{0,uj;m) 

We first estimate the imaginary part. From Eq. (C5), we have 

Im^(0,.;™).,o..^ (-/(I) + .(!)) . Ini.). (C7) 
Here we have utilized the inequality w 3> m. Using the approximate formula 

n{x) ~ T/x for a: < T, (C8) 

we get 

ImB(0,L^;m)T^o-^-. (C9) 

The real part is given by Eq. (C3). To evaluate it, it is found convenient to introduce an intermediate scale T* 
which satisfies the relation T ^ T* ^ a; so that T* separates the hard region from the soft region. The contribution 

from Ef > T*, which is denoted by KcB{0,lu; 'm)^[^'^ , is evaluated to be 

* 1 1"°° 1 

ReB(0,c^;m)^^>^ ^ P j^^ dEf— {f{Ef) - n{Ef)) . (CIO) 

Here we have used the inequalities lo ^ m and Ef > T* ^ u in the hard region. We get the following expression by 
integrating by part: 

• 1 / roo ^ ft \ 

ReS(0,c.;m)^^f ^ - — ^- [f [E f) - n{E f)]'^, + V dEf — —Ef{f{Ef)-n{Ef))j. (Cll) 

By inserting Eq. (C8) and f{x) ~ 1/2 for x ^ 1, we see that the first term is of the order of T/T*. It is easily 
confirmed that the second term is of the order of ln(T*/r) by performing partial integration again, so the second 
term can be neglected; 

ReB(0,a.;m)^g^*c.^|l. (C12) 
On the other hand, we have for the contribution from Ef < T*, which is denoted by Re^(0, cj; 'm)^^_^^ : 

Here again, Eq. (C8) has been inserted and the inequality uj ^ m has been taken into account. The term which has 
f{Ef) has been neglected since it is much smaller than n{E\,). Performing partial fraction decomposition and utilizing 
the inequality w 3> m, we get 

Rei?(0,.;Har--^ + ^- (C14) 

We see that the imaginary part, Eq. (C9), is much larger than any term of the real part, Eqs. (C12) and (C14). 

A remark is in order here: In the present analysis, we have considered only the leading term in the high temperature 
expansion as in [14, 17], where the simple approximate formula n{x) ~ T/x is adopted. In fact, all the analyses given 
in this Appendix will be based on this formula. Of course, one could include the next-to-leading term by using the 
formula n[x) ~ T/x — 1/2, which would lead to logarithmic terms derived in [20, 21]. The inclusion of such terms is 
beyond the scope of this work. 
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b. B°{0,Lj-m) 

The imaginary part of i3°(0, w; to)t/o is given by Eq. (C6), which is estimated to be 

ImSO(0,L^;m)T^o- T^. (C15) 

The derivation of this expression is the same as that of Eq. (C9). 

The real part is given by Eq. (C4). We introduce T* which satisfies T ^ T* ^ lu, as before. The contribution 

from the Ef > T* region, which is written as Rei?°(0, cj; m)^^p^ , is evaluated as follows by inserting Eb ~ Ej and 
Ef ^ uj ^ m again: 



ReB''(0,c.;m)^^f J^^ dEfEf-^{f{Ef)+n{Ef)) 

--^[[ d^f - [ ^^/) ^/(/(^/) + -(^/))- 

Inserting Eq. (C8), we obtain 

7^2 1 j-T' 

RcB\Q,uj-mp>^ ^~ — + —^ cLEfT 
^ 16a; A-K^LJ J„ ' 



(C16) 



(C17) 



16a; Att'^uj 



The contribution from Ef < T* region, represented by Re^*'(0, oj; 771)^^0''" , is estimated as follows by using 
Eq. (C8) and the inequality lu ^ m: 

Utilizing the inequalities T ^ T* ^ lu, Rei3"(0, oj; m)^^^"^ is evaluated as follows by utilizing partial fraction 
decomposition: 

ReBO(0,a;;™)a<- - ^ (— + ' (^l^) 

Since the largest term is — r^/(16a;), we have 

RcB°(0,a.;TO)T^o--^, (C20) 
loo; 

which is much larger than the imaginary part. A few remarks are in order here: the last formula leads to the well- 
known result of the HTL approximation in QED. It is easy to understand if we remember that we can neglect m and 
apply the HTL approximation in this situation, T ^ uj ^ m. 

From Eqs. (4.12) and (4.13), it seems that Tm/ [Attuj) in Eq. (C19) yields terms which are of the order of g^Tuj/m 
and g^Tuj^/a/m in E+(0,a;), which implies that the massless limit (ra/T — > 0) can not be taken. However, from 
Eq. (C14), we see that these terms are canceled by terms which come from Eqs. (4.8) and (4.9). 

The results of this section, Eqs. (C12), (C14), (C17), and (C19), reproduce leading term of asymptotic form of 
some integrals used in [21] since the QED limit, m/T, m/o; — 0, is included in the limit considered in this section, 
T » oj ^ m. Results in other sections are not related to them because their situations are far from QED. 
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2. T > m > 

a. B{0,ui;m) 

The imaginary part of 5(0, a;; rn)T^o is given by Eq. (C5), which is estimated as fohows: 

f 'W?\ ( rn? \ \ — TO^ / m? 
— + n — = ^ exD 



iJS^cxpf^-^j (m2»^Tcase) 
-gij (m^^wTcase) 



(C21) 



In the first hne the inequahty m ^ uj has been taken into account, while in the second hne, either 3> ujT or 
<C wT has been used. 
Next, we estimate the real part by separating the cases ^ uiT and ujT S> m^. 

• ^ case 

We begin with Eq. (C3). Wc again introduce an intermediate scale T* satisfying T ^ T* ^ m. The contribution 
from Ef > T* is evaluated as follows: 

• 1 

Re^(0,w;m)^^>^ - j^^ dEfEf{S{Ef) + n{Ef)), (C22) 

where the relation Et,:^Ef^m^uj has been utilized and the contribution from Ef > rn^/uj has been neglected 
because it is suppressed by the Boltzmann factor ^ exp{—m'^/{ujT)). Dividing the interval of the integral as in 
Eq. (C16), the following expression is obtained: 

ReB(0,w;m)^^>^* ^ - If dEj - [ dEf] Efif{Ef) + n{Ef)) 

^TT m Jo J (^23) 



In the last line, we have inserted Eq. (C8). On the other hand, using Eq. (C8) and m ^ a;, the contribution from 
Ef < T* is estimated to be 



Utilizing the inequalities T* <C m /u and m ^ lo, we obtain 

ReB(0,c^;m)^^f' ^ - + (C25) 

which is obtained by partial fraction decomposition. Taking only the largest term, we have 



2 



ReB{0,iu;m)T^o--^, (C26) 



which dominates the imaginary part. 
• TO^ ^ ujT case 



contribution from Ef > T 



We introduce T*, which is an intermediate scale satisfying T ^ T* ^ m?/u!. From Eq. (C3), we have for the 

RcB(0,c.;m)^^f c,^-^ j^^ dEf — {f{Ef)+n{Ef)) 



(C27) 
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In the first line, we have used the inequalities m ^ uj and c:^ Ef > T* ^ m^/w, while in the second line an 
approximation similar to that after Eq. (Cll) has been utilized. 

Substituting Eq. (C8), the contribution from Ef < T* is estimated to be 

,Ef<T' m^T Ej 1 



Tfl T f 

ReB(0,c^;m)^^<^* ~ - dEj 

m^T T 



Ej + m^m^-AEj.^ (C28) 



In the second line the inequality T* 3> m? /ui was used. Now we see that every term in the real part is much smaller 
than the imaginary part. Combining these results, we have 

\ {—hi — Y (m^ > wT case) 

B{^.-:^)T^o^[jh^ K«..Tcase . (^^9) 

b. B"{0,uj;m) 

First, let us estimate the imaginary part of i?°(0, ui; rnjT^o given by Eq. (C6): 



~n, . -I " j cxp - (m^ > case) 



(m2<a;rcase) . 



We have arrived at this expression in the same way as in Eq. (C21). 
• ni^ ^ ujT case 

In this case, the real part of i3°(p, w; m)T=io which is given by Eq. (C4) is evaluated as follows. Using an intermediate 
scale T* satisfying T ^ T* m, we have for the contribution from Ef > T* 

roo 

ReB°{0,Lo;mp//^^-^ dEfE){f{Ef)+n{Ef)), (C31) 

TT Tfl J J^* 

where the contribution from Ef > rn^/uj has been neglected since it is suppressed by the Boltzmann factor ^ 
exp(— TO^/(a;T)). Also the inequality m 3> w has been taken into account. Dividing the interval of the integral as in 
Eq. (C16), we obtain 

ReB°(0,a.;m)^g^* ( H dEf - f dEf] E}{f{Ef) + n{Ef)) 



*3 



8 m4 37r2m4' 

In the last line, we have substituted Eq. (C8). 

Next, we estimate the contribution from Ef < T* . Since the leading terms cancel out, we retain in Eq. (C4). If 
we define ilj. = ± m^, the contribution from Ef < T* is estimated as follows: 



Rei^0(O,.;™)-^r.gp^ dEf^r^,^^^^^,, (C33) 



where we have used Eq. (C8). Performing partial fraction decomposition, we have 

ReB° 0,a;;m / dEf 1 ^ -—^ ^ -j— — — 

^' 47r2c^ h ^\ E]+m^nl Au^Cii Ej - nt/{4uj^) J ^^^^^ 
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In the second line the inequahty T* ^ m has been utihzed. Notice that because of cancehng out of the leading terms, 
we have to retain the next-to-leading term in the expansion of ln(l -I- 2wT*/ri^). The largest term is the first term 
of Eq. (C32), so 

Re^O(0,w;m)T^o - (C35) 

which is found to be much larger than the imaginary part. 
• 771^ ^ uT case 

The real part is given by Eq. (C4). As is before, we introduce an intermediate scale T* satisfying T ^ T* ^ 
rr? jio ^ TO. Then 

R,BO(0.„;Hftf.-ll + g:. (C36) 

This expression has been obtained in the same procedure as in Eqs. (C16) and (C17). 

The contribution from < T* is estimated as follows. Using Eq. (C8) and the inequality to 3> i^, we have 

ReB°(0,a;;ra)f;^V ^^^P / dEfT „ ^ , — ^ 

^ 27r2 io ^ £:| + to2 ,^4 _ 4^2^2 ^^^^^ 

TT* m^T 



In the last line we have utilized the inequality T* ^ to^/w ^ ?7i. Since T* S> / [uP'T) , the largest term is found to 
be — T2/(16a;), so we have 

ReB°(0,w;TO)T^o ^ (C38) 

This is much larger than the imaginary part. 
Combining these results, we have 

fSO/r. ^ I (to2 > case) ,r^^r^-. 

B"(0,^;TO)r^o^ <^ 2^ ^ <^39 

I -fb (m2<cjrcase) . 



3. m > T > oi 
a. B{0,uj\m) 

We start with an estimate of the imaginary part of B{0,lo; m)T^o given by Eq. (C5). For to ^ T ^ w, we have 

77^2 / y^2 \ 

ImB(0,.;TO).,„.-^exp(^-_j. (C40) 

The derivation of this expression is almost the same as that of Eq. (C21) in the to2 ^ uiT case. The real part of 
B{0,uj;'m,)T^o given by Eq. (C3) is estimated as follows: 

1 r°° -to2 
ReB(0,.;TO).^o - ^P [ dEfEffiEf) ^,^^^,^, . (C41) 

Here we have used the inequality m uj and n{Eb) has been neglected because of the suppression factor by 
exp(— to/T). Notice that the contribution from Ef > m^/ui 3> T is negligible owing to the suppression factor 
by f{Ef). Then, we see that the denominator of the integrand can be approximated as to^ — AE'joj'^ ~ to^. Thus we 
have 

ReB(0,^;TO)T^o ^ ^^P dEfEff{Ef) = (C42) 
which is found to dominate the imaginary part. 
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b. B°{0,Lj-m) 

We estimate the imaginary part of B'^{0,uj; m)T/Oi which is given by Eq. (C6). Utihzing the same procedure as in 
Eq. (C40), we have 

Imi^"(0,.;™),,o-^cxp(-£;). (C43) 

The real part, which is given by Eq. (C4), is estimated as fohows: 

RcS"(0,c.;m)T^o ^ f Ej f {E f) -^-^^ . (C44) 

Here we have utiHzed the inequahty m ^ lo and neglected n{Eb) since its contribution is suppressed by the factor 
exp(— m/T). Following the procedure used in obtaining Eq. (C42), we have 

ReBO^O,.;m)r,o / dE,E^,fiE,) = (C45) 

One sees that this term dominates lmi?'^(0, a;; to)t5^o- 

4. T > oj ~ m 

a. B(0,u};m) 

We start with an estimate of the imaginary part of _B(0, a;; to)t^Oj which is given by Eq. (C5). Utilizing the 
inequality T ^ uj ^ m and substituting Eq. (C8), we have 

lmi?(0, oj- m)T^o ^ ^ J— ^ - -■ (C46) 

In the evaluation of the real part, wc introduce T* satisfying T 3> T* 3> w. From Eq. (C3), we arrive at the 
following expression by using the inequalities Ei, ^ Ef uj, m: 

RcB{0,^;m)pJ'/ -^^P j^^ dEf— {nl f {E f) - n\niE f)) (C47) 

with = cj^ ± m?. In the last line we have performed the partial integration and the order estimate as we did for 
Eq. (Cll). 

Inserting Eq. (C8), the contribution from Ef < T* is estimated to be 

Using the inequality T ^ T* ^ u ^ m and carrying out partial fraction decomposition, the following estimation is 
obtained: 

We see that both of the real part and the imaginary part are of the order of T /uj\ 

B(0,c^;to)t^o- -• (C50) 
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B°{0,ij;m) 



We make an order estimate of the imaginary part of i3''(0, w; to)t/o given by Eq. (C6). Repeating the same 
derivation leading to Eq. (C46), we have 

ImBO(0,a.;m)T^o =^ J /"''^TT ^ ^- (^51) 
From Eq. (C4), we estimate the contribution from the Ef > T* region to its real part as follows: 

ReB»(0,.;,„)?S-.-|l + |^. (C52) 

Here we have introduced T* which satisfies T ^ T* ^ m, and used the same approximations as that adopted in 
Eqs. (Cf6) and (Cf7). 

On the other hand, the contribution from Ef < T* is estimated as follows by using Eq. (C8): 



ReS°(0,w;m)^lo^* ~-^P / dEfT 



E^f + m2 f7i - 4£;2w2 
/ X (C53) 

T f TT fnl \ f 1\ ^2 2^2 / 1 ^ ^ ' 

T* -m- + — ^ - 1 + m^n^ — i ^ m- 



47r2cj \ 2 2 \2uj^ J\ T* J ~ nt + im'^co'^ \ 2 T 

with = cl|2 — m^. In the last line, the inequality T* ^ uj ^ ra has been utilized and partial fraction decomposition 
has been performed. Since the largest term is —T'^/{1Qlu) in Eq. (C52), we have 

ReB°(0,a;;m)T^o--^, (C54) 
low 

which is much larger than the imaginary part. 
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